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Abstract 

m 

Background: It has been observed in heavy ion collisions at both the Relativistic Heavy Ion 
Collider (RHIC) and the Large Hadron Collider (LHC) that the high energy jets are strongly 
quenched because of its interaction with the created hot colored medium. Yet the azimuthal 
anisotropy of the quenching has not been studied thoroughly so far, which should provide 



important information about jet-medium interaction and the medium's properties. 



Purpose: We present with full details a systematic quantification of the hard probe azimuthal 
response to the geometry and fluctuations of the hot QCD matter at both RHIC and LHC. 
We also examine the hard-soft di-hadron correlation arising from the hard and soft sectors' 



C\| 1 responses to the common fluctuating initial condition. 



Methods: An even-by-event MC simulation is employed. Different geometrical jet-quenching 
models are tested. The azimuthal anisotropy of jet quenching is extracted and decomposed 
as harmonic responses (for n=l-6) to the corresponding harmonics in the initial state. 



Results: We show that such jet response harmonics are sensitive to the jet quenching models 
^ 1 as well as to the bulk matter initial compositions. Their centrality dependence puts a 

strong constraint on the path-length and medium-density dependence of jet energy loss. 
The computed hard-soft di-hadron correlation shows a strong peak on the near side in RHIC 
central collision. The triggered correlation in the non-central collision is also presented. 

Conclusions: Only the jet-quenching model with near-Tc enhancement survives the 2nd- 
harmonic-test by RHIC and LHC. Other harmonics in this model are consistent with the 
available data. We also demonstrate that the experimentally observed "hard-ridge" can be 
explained in our calculation and that its trigger-azimuthal-angle and associate-^ dependence 
could also be qualitatively understood. 
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I. INTRODUCTION 



In this paper we present with full details a systematic quantification, on an even-by-event 
basis, of the hard probe response to the geometry and fluctuations of the hot QCD matter 
created in heavy ion collisions at both the Relativistic Heavy Ion Collider (RHIC) and the 
Large Hadron Collider (LHC), substantiating previous studies reported in [lH^]. In the high 
energy Nucleus-Nucleus (AA) collisions, the hot deconfined QCD matter, so-called quark- 
gluon plasma (QGP), is believed to be created. Also created in such collisions are jets from 
the initial hard scatterings: such a hard probe penetrates the soft matter and eventually 
converts into hadrons. Along its path through the medium the jet keeps interacting with 
the medium and losses energy. The quenching of jets due to the energy loss can be used to 
characterize the jet-medium interaction as well as the medium properties, which provides 
an imaging tool of QGP often called jet tomography (see reviews in Refs. [5|). 

In this study we focus on the single hadron observables, in particular for the high p t region 
dominated by jet physics. The traditional observable for quantifying the jet quenching is 
the so-called nuclear modification factor Raa, i-e. the ratio between the hadron production 
in AA collision and the production in NN collision scaled by the expected binary collision 
number. The Raa has been measured to be substantially less than 1 for central collisions 
first at RHIC and now also at LHC, indicating a strong overall jet quenching or equivalently 
that the hot medium is rather opaque to colored hard probes. While the hot medium is 
opaque, the interesting question now is whether the medium at LHC is more or less opaque 
as compared with RHIC (apart from the trivial density increase). It is interesting to notice 
that models with one parameter characterizing the medium opaqueness (or jet-medium 
interaction strength) and tuned to describe RHIC central Raa seem to consistently show 
over suppression as compared with LHC data and a sizable reduction of the opaqueness from 
RHIC to LHC is strongly indicated @, @, 0. 

The central issue to be addressed by this study is the azimuthal anisotropy of the jet 
quenching, i.e. the dependence of Raa on the azimuthal angle <fi (relative to the event 
plane angle ijj EP ) 0, @-15|. Such anisotropy originates from the asymmetric spatial distri- 



butions both of the bulk matter and of the initial hard collisions due to both geometry and 
fluctuations. 

The anisotropy from geometry has been known and studied for long time. In most AA 
collisions, the overlap zone has an almond-like geometry in the transverse plane, which is 
then inherited by the created medium at early time. As the result, the jets with different 
transverse orientations (at mid-rapidity) would see different medium "thickness" and lose 
different amount of energy, which leads to the azimuthal anisotropy of final-state high p t 
hadrons. The dominant component of this anisotropy lies in the 2nd order harmonic (often 
called the elliptic component) in Raa(4>) ~ 1 + v\ cos[2(</> — ip EP )]. This v% should be 
differentiated from the commonly known elliptical flow i>fj because of their distinctive origins: 
one is from jet quenching and the latter is from medium's collective expansion. At RHIC, 
most model calculations predicted a v\ much smaller than what the data indicates. To 
resolve this, a nonperturbative model was suggested in Ref. |l( (named as NTcE model later) 
with a more radical insight that the jet-medium interaction bears non-trivial dependence 
on matter density and is strongly enhanced in the near-T c matter, as motivated by the 
"magnetic scenario" for sQGP [16l. Il7|. As we will discuss later, this model naturally expects 
a reduction of the jet-medium interaction from RHIC to LHC due to the nearer-to-T c nature 
of RHIC matter and a shift toward higher density matter at LHC. 
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One significant progress over the past two or three years, came from the realization \lt 
and subsequent detailed investigation of strong fluctuations in the initial condition of heavy 
ion collisions. The above almond-shape picture for v\ and v\ assumes a smooth geometry for 
the medium. The reality turns out much more complicated and interesting: event by event, 
the nucleon positions in the two colliding nuclei fluctuate, leading to the so-called initial 
state fluctuation (ISF), i.e. the fluctuation of the medium's initial (entropy and energy) 
density distributions, which are then translated by hydrodynamic expansion into different 
harmonic flows in the finally observed soft hadron production. These have been observed 
in the RHIC and LHC experiments 19M23j and also demonstrated in various hydrodynamic 



calculations [18|, |24J-|371]. In addition, the di-hadron correlation seems also related: due to 
the existence of other harmonics v s n in the collective flow besides the elliptical flow, the 
di-hadron azimuthal correlations get a contribution from simply folding together the two 
single particle distributions and showing a strong peak on the near side ("soft ridge") as 
well as nontrivial structure on the away side depending on the centrality and hadron p t (e.g. 

SSHHSisI). 

Then how about the jets? Following the same argument about induced by the almond- 
like geometry, the high p t azimuthal anisotropy can be affected also by the initial anisotropy 
from strong event-by-event fluctuations in the geometric distributions of both the soft matter 
density and the initial jet spot profile (the collision density). Some earlier discussions and 
data from RHIC and LHC can be found in Refs. 13, BEE). In Ref. Q we've made the 
first attempt to quantify systematically the high p t harmonics as hard probe response to 
fluctuating initial conditions, for RHIC central collisions with the fluctuations implemented 
via cumulant expansions. Here we continue and significantly advance the study along this 
new direction about the connection between ISF and jet azimuthal anisotropy (z/„) as initi- 
ated in Ref. [3|, by fully implementing the event-by-event Monte Carlo initial conditions and 
extending the quantification to all centralities for both RHIC and LHC collision energies. 
We emphasize that this is an important shift for jet quenching physics, i.e. from studying an 
average anisotropy pattern Raa(4>) ~ 1 + ^2 cos[2(0— ^ EP )] toward event-by-event extraction 
of full anisotropy information Raa^) ~ 1 + J2 n =i2 3 ^cos[n(0 — ip^)]- m particular, as 
alluded above, we want to gain a coherent picture about the jet quenching and its azimuthal 
anisotropy from RHIC to LHC, as already partly investigated in Ref. [4(. The hope here 
is that by using geometric data at both collision energies, a strong constraint can be put 
on different jet quenching models, especially the energy loss dependence on the path length 
as well as the medium density (temperature). In addition, understanding the jet response 
anisotropy provides a complimentary and sensitive tool beyond the bulk hydro response for 
probing the initial conditions. The analysis in the present paper will further address a num- 
ber of important questions, including: the event-by-event determination of eccentricities e n 
and participant angles ip n at both RHIC and LHC experiments; the jet azimuthal anisotropy 
in terms of harmonics arising from geometry and fluctuations; comparison of the jet har- 
monic responses in different energy loss models; the transfer of eccentricities to jet response 
haromincs and the angular correlations between ip^ and ip n ; examination of the separate 
contribution to the anisotropy from initial jet spot profile fluctuation and from the medium 
density profile fluctuation; the sensitivity of the jet response anisotropy to the composition 
of the matter density (in terms of participant and collision profiles); etc. 

Finally, in parallel to the study of di-hadron correlation in the kinematic regime domi- 
nated by soft hadrons, we will also study the azimuthal correlation between hard and soft 
hadrons [IB], |46[ (hard-soft correlation) . In particular we examine the contribution from their 
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common correlations to the fluctuating initial conditions, as also initiated in Ref. j3|. Our 
improved calculation confirms what we previously found for the central collision at RHIC: 
the hard-soft correlation shows a strong peak at near side ( "hard- ridge" ) and double- hump 
structure on the away side due to the concurrent azimuthal harmonics in both soft and hard 



particle distributions, which is a possible explanation of the "hard- ridge" [47H50j . Recently, 



STAR collaboration has extracted the trigger-azimuthal-angle dependence of the di-hadron 
correlation 51-53 . In Ref. the trigger-angle (relative to -?/> EP ) and p t dependence of all 
the harmonic components extracted from the di-hadron correlation was studied from the hy- 
drodynamic calculations which should be good in the kinematic regime of associate hadron 
Pt around 0.15 ~ 3 GeV and trigger p t around 3 ~ 4 GeV). However for even higher p t 
trigger (e.g. 4 ~ 6 GeV), the picture should be different as jets become the dominant source 
and the correlation is of the hard-soft type, which we will study based on our computation 
of jet response to ISF. We will also re-derive the formula for the triggered correlation [|34j 
by lifting the often-used (and not fully correct) assumption that different event plane angle 
ip£ (for hard hadron) and ip^ (for soft hadron) are totally correlated with ip EP (i.e. ?/>f). 

Before getting into the main body of this study, we would like to emphasize two important 
discussions on the near-T c enhancement model of jet quenching that have been included as 
the Appendix C and D. The first issue is about the possible relation between the strong 
near-T c peak seen in the QCD trace anomaly (a measure of non-conformal behavior) and the 
proposed near-T c enhancement that arose during a discussion j55|, and here we will provide 
an answer in Appendix C to show that indeed the two phenomena could be simultaneously 
and consistently understood from the scenario of a near-T c plasma of magnetic monopoles. 
The second issue concerns a quantitative estimate of the evolution of medium opaqueness 
from RHIC to LHC: while it was pointed out before that the strong near-T c enhancement 
of jet-medium interaction naturally expects a less opaque medium at LHC, in the Appendix 
D it is quantified to be a reduction of about 30% for the average jet-medium interaction, 
i.e. < k >rhic:< k >lhc~ 1 : 0.72 — in consistency with a number of recent jet quenching 
analyses |6|, |7J, l56l - 58j. These are two important points and we choose to leave them in 



the Appendices primarily because they are somewhat less attached to the stream of main 
discussions in this paper. 

The rest of the article is organized as follows. In Sec. [TTJ we collect the details of our MC 
simulation. The initial state is modeled by using Glauber model. Three different geometric 
jet energy loss models are presented. The structure of the simulation code can be found 
in Appendix |A] Sec. IIHI summarizes our results for the jet anisotropy at RHIC and LHC 
and the comparison of different models. In Sec. HVl we discuss the (un)triggered hard-soft 
correlations. The re-derived formula can be found in Appendix [B] A summary and some 
discussions will be given in Sec. |V] 



II. MC SIMULATION 

In this section, we discuss the details of the MC simulation. The structure of the code 
can be found in Appendix [A] 
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Po (fm 3 ) 


Ra (fm) 


d A (fm) 


sfB (fm) 


cj inc i (fm 2 ) 


So (fm" 3 ) 


n 


RHIC (0.2 TeV) 


0.170 


6.42 


0.45 


0.544 


4.2 


116 


8.1 


LHC (2.76 TeV) 


0.166 


6.67 


0.44 


0.660 


6.2 


291 


6.0 


LHC (5.50 TeV) 


0.166 


6.67 


0.44 


0.680 


6.6 


364 


5.4 



TABLE I: Parameters in the simulations for both RHIC and LHC. The explanation and the refer- 
ences for them can be found in the text. 



A. Glauber model 

In the AA collisions, event by event, the positions of the nucleons fluctuate, as the result 
of the "measurement" of their positions (in the quantum language). This leads to the 
fluctuation of matter density, i.e. ISF. In the simulation, it is realized in a simplified way: 
sampling every nucleon position in the nuclei according to its density as measured in the 
low-energy scattering process [HH . The short range correlation between nucleons is included 
via requi ring a smallest distance between them, which is set 0.4 fm in most simulations 



30j, 1591 ) including this work. The Glauber model is applied to deal with multiple AW 



collisions [59] : binary collision happens only if the transverse distance between two nucleons 
is smaller than a/ cr- mc i/iT. Here o"i ne i is the total inelastic AW scattering cross section at 
the corresponding center mass energy. The trajectory of each nucleon is always along the 
beam direction. The parameters in the simulation of the initial state are summarized in 
the following. First, the density in each colliding nuclei is parameterized in terms of the 
Wood-Saxon form [60j: 

ex p[ —if 1 ] + 1 

The values of oq, Ra, and (Ia for Au and Pb nucleus as used in RHIC and LHC can be found 
in Tab. II l - 

For the density in a nucleon, a Gaussian distribution is assumed in the 



transverse plane [60 



/ \r A 

e M-— B 



TJr ) = ^ 2B ' , (2) 
pV ; 2ttB ' V ; 

where B value can be found in Tab. |T] for three different energies. Second, the total cross 
sections cr inel for AW collisions at these energies are also listed in Tab. [B After the binary 
collisions are sampled by using the Glauber model based on a inc i, we sum up contributions 
of all the "wounded" nucleons to get the participant density in the transverse plane: 

P P (r ± ) = J2 T P (\r ± -ri\) . (3) 

i=i 

Here N p is the number of wounded nucleons. Similar procedure can be applied to calculate 
the collision density p c (r ± ). 



The finite size of nucleons is taken into account when fixing the parameters in p(r) as the sampling 
probability. 
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To make a connection between the initial state and the state of equilibrated medium, we 
assume that the entropy density S at r = tq = 0.6 fm/c (around the equilibration time) is 
proportional to p p at RHIC (0.2 TeV) [IrS], and a two-component profile, (1 — 5) x p p /2 + 5 x 
p c , at LHC (5 = 0.118 for both 2.76 and 5.5 TeV) (61 6^] . The maximum entropy density 
in the central collision, Sq, is listed in Tab. [I] for different experiments [62-64j]. which are 



calibrated to the observed multiplicities in the hydrodynamic calculation. Based on this, 
we compute the proportionality coefficient in the relation between p p and p c , and S(r ). 
Before To, we assume effectively for jet quenching, the entropy density S grows linearly with 
time, while after that, the longitudinal expansion is applied and S decreases as 1/r. See 
Ref. (39| for a detailed discussion about the effects of different pre-equilibrium models on 
jet quenching. Unfortunately, the medium's transverse dynamics is simplified without its 
hydrodynamical evolution, i.e. the shape is frozen, which will be discussed later. For the 
jets, their productions at the early stage of the collision should distribute according to p c 
and are isotropic in their momentum space. In the following section, we focus on the jet 
energy loss when traveling through the medium. 



B. Jet quenching models 

Three different geometric models are applied to compute the jet energy loss, which have 
distinctive geometric features (e.g. the path-length dependence) and matter-density depen- 
dence that are most crucial for describing geometric data [lTM fill 12, 39, Suppose the 



jet is produced with initial energy Ei. After traveling a path P in the medium, the ratio 
between its final energy Ef and the initial energy E iy i.e. the suppression factor fp, is given 
by 

/ P = exp(- f K[s(l)}s(l)l m dl) . (4) 



In this expression the s(l) is the local entropy density at a given point on the jet path, and 
the k(s) is the local jet quenching strength which as a property of matter should in principle 
depend on the local density s(l). There can be different choices of the parameter m for 
path- length dependence (e.g. LPM- motivated quadratic or AdS/CFT-motivated cubic) and 
of the jet-medium interaction k(s). In this study, we compare the following models the 
near-T c enhancement (NTcE) model, L 2 model, and L 3 model. The NTcE model assumes 
m — 1 and introduces a strong jet quenching component in the vicinity of T c (with density 
s c and span of s w ) via 

k(s) = «q[1 + f exp(-(s - s c ) 2 /s 2 J] , (5) 



with £ = 6, s c = 7/fm 3 , and s w = 2/fm 3 . (see [lj for details.) In the L 2 (L 3 ) model, 
k(s) = Ko is constant and m — 1 (m = 2). In each model, the parameter k,q will be fixed by 
Raa « 0.18 in the - 5% collisions at RHIC y/s = 0.2 TeV, and then be applied at LHC 
= 2.76, 5.5 TeV. 

To calculate the azimuthal angle dependence of the nuclear modification factor, we can 
apply the following formula: 

Raa(<P) =< (fp) n - 2 >p W ■ (6) 

Here < >pu) means averaging over all jet paths with the propagation orientation fixed at 
angle (ft relative to the reaction plane, including all the possible initial jet production points 
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[distributed according to p c (r J -)]. The exponent n comes from measured high p t spectrum 
in reference p-p collision, which depends on collision energy and is shown in Tab. [J for the 
three different energies (65h67l|. To calculate the over all modification factor, Raa, we need 
to average over all the jet paths with different orientation angles that are equally distributed. 

Moreover, in these models, assuming the fractional energy loss leads to the Raa{^>) in- 
dependent of p t . This at RHIC energy may be justified by the approximate "flatness" seen 
in the Raa vs. p t data. At LHC energies, we expect the p t dependence of the azimuthal 
anisotropy, i.e. Raa(4>) I Raa, to be weak, which should be dictated by the length and density 
dependence of the jet energy loss in the transverse plane. The results from such modeling 
apply only to the high-pt region, e.g. p t > 6 GeV at RHIC 0.2 TeV, and p t > 8 GeV at LHC 
2.76 TeV 2 . 

In the simulation of each event, we apply the Fourier-decomposition to the Raa{4") I Raa- 



Raa(4>) = Raa 1 + 2 5>n cos[n(0 - ^)] 



(7) 



In each event the azimuthal anisotropy can be represented by the collection of and i\) J n . 
The importance of different harmonics decreases with increasing n except the 2nd harmonics, 
as will be shown later. So we focus on the first six harmonics in this work. 



III. THE ISF AND JET AZIMUTHAL ANISOTROPY 



A. The ISF 



To quantify the ISF in each AA collision, the following eccentricities e n are defined, as 
motivated by the cumulant expansion method in Ref. 29| : 



< r n cos[n(0 — ipn)} > 



<r n > 



arctan 



n 



< r sin [nc 



> 



< r n coslnd)) > 



71 

n 



< r 3 cos[(0 — ^i)] > 



< r 6 > 



arctan 



< r 3 sin(0) > 

< r 3 cos(0) > 



7T 



(8) 
(9) 
(10) 
(11) 



The ip n is the reaction plane angle of nth harmonics. < • • • > means averaging over the 
entropy density at r . 

Fig. [1] shows the e n spectrum for b = fm (most central collision) and b = 7 fm (peripheral 
collision) at three different collision energies. For all the experiments, in the central collision, 
e n are close while in the peripheral collision, dominates due to the almond like geometry. 
Moreover, e n at LHC 2.76 TeV and LHC 5.5 TeV are almost the same. In Figs. 121 and fT5| 
we plot the distributions of ipn relative to the reaction plane for 6 = and 7 fm in three 
cases [29j]: in the central collision, ip n do not have specific orientation; in the peripheral 



2 At LHC 5.5 TeV, the threshold maybe somewhat higher. 
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FIG. 1: (Color online). The spectrum of e n at b = 0, 7 fm for RHIC with y/s = 200 GeV and for 
LHC with y/s = 2760, 5500 GeV. 
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FIG. 2: (Color online). The distribution of ip n in b = 0, 7 fm for both RHIC (y/s 
Similar results for LHC (y/s = 2760, 5500 GeV) can be found in Fig. M 



200 GeV). 



collision, tpi and ipi tends to lie in the 7r/2 and ir/A directions; ip2 is strongly correlated 
with the reaction plane (in the plot, ?/Vs probability density is scaled down by a factor 0.3); 
other odd ip n distribute randomly. We also check that the correlation between two angles 
vanishes in the central collisions. Figs. |3] and [16] list the b dependence of the eccentricities, 
which shows the similarity in three experiments and also the emergence of the 2nd harmonic 
dominance in the peripheral collisions. In addition, Fig. H] shows the relations between N p 
(N c ) and b. The error bars indicates the r.m.s of them, which are small except in the very 
peripheral collisions. In the following, we will show the results against b. Fig. H] can be used 
to compute the corresponding N p dependence. 



B. The jet azimuthal anisotropy 

In this section, we focus on results for the jet azimuthal anisotropy. Figs. 151 and [171 show 
the 1;% spectrum [see Eqn. ([7])] in both most central collision (6 = fm) and peripheral 
collision (6 = 7 fm) as calculated by using different jet energy loss models. We can see that 
the 2nd harmonics dominates in the peripheral collision reflecting the almond-like geometry 
in the peripheral collision, and other harmonics decrease with increasing n in both central 
and peripheral collision. More importantly, for the first three harmonics, the NTcE and 
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FIG. 3: (Color online). The b dependence of e n for both RHIC (y/s = 200 GeV). Similar results 
for LHC (y/s = 2760, 5500 GeV) can be found in Fig. M 
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FIG. 4: (Color online). The b dependence of N p and N c for both RHIC (y/s = 200 GeV) and LHC 
(y/s = 2760, 5500 GeV). 



L 3 result for RHIC (0.2 TeV) are close, while for LHC (2.76, 5.5 TeV) the NTcE and L 2 
result are close. The b dependence of u% in the three experiments as predicted by the 
models is shown in Figs. El [7J and [181 (curves labeled as "model i/"). < cos[n(^ — ip n )] >, 



which reflects the correlation between i\) n and ip n 68] , is also plotted in these figures against 



different h for different harmonics (curves labeled as "model c"). As illustrations, Figs. dHl 
I2"0"l and HU attached in the end of this article, show the distribution of ipn — i>n computed 
by using different models for 6 = and 7 fm in the three experiments. Interestingly, the 
authors in Ref. [30] presented similar distribution of — ip n [Here ip^ is the angle of nth 
harmonic event plane reconstructed from the final low p t hadron azimuthal distribution; see 
Eq. fIBlj) ]. By comparing the two, we find that for the 2nd and 3rd harmonics, ip^ and i/j^ 
have similar distribution relative to the ip n , although ip^ distributions computed by three 
models differ somewhat; for the other higher harmonics, an d ipn gains much broader 
distribution relative to ip n , and < cos[n{ip^ — ip n )} > changes from positive to negative with 
b decreasing (the value of the transition b depends on individual jet energy loss model). 
On the other hand, the experimental measurement of u% is normally not the u% shown in 
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FIG. 5: (Color online). The spectrum of for RHIC (y/s = 200 GeV) and LHC (y/s = 2760 GeV) 
at b = 0, 7 fm, based on three different models. Similar results for LHC (y/s = 5500 GeV) can be 
found in Fig. [T71 



Figs. El and[T8| but its projection to ip%, i.e. < ufc cos[n(i/r[ — ip^ )] > 3 . Because of missing 
soft dynamics in our simulation, we can not calculate < cos[n(/0^ — ip^)] > a ^ ^ ms stage. 
However, for the 2nd harmonics u!j\ based on the < cos[2(/02 — ^2)] > shown here and the 
■02 — ^2 distribution shown in Ref. [30|, the effect of this angle dispersion should be less than 
5% reduction for peripheral collision (6^6 fm). This provides solid ground for comparing 
the data for v\ with our results to differentiate different jet energy loss models, as carried 
out in Q . There it was found that the RHIC v\ data favor the NTcE and L 3 models, while 
the LHC (2.76 TeV) v\ data favor the NTcE and L 2 models. This is consistent with v\ 
shown in Figs. El [7J and [18j Meanwhile, v\ and v\ delivers valuable information about the 
order of the magnitudes of the measured harmonics, although the angle dispersion can bring 
substantial reduction yet without changing the sign. It is conceivable that the experimental 
data on will bring important constraints to the jet quenching models. A tentaive try 
along this direction is pursued in Ref. Q ■ 

Let's also take a look at the correlation between ip^ and the participant plane ip2- For 
the central collision, we can expect ty n to distribute randomly w.r.t ip2, which is confirmed 
by our simulation results for RHIC (0.2 TeV). However, in the peripheral collisions, things 
become tricky. Fig. [8] shows the ip n — ip2 and ip£ — ip2 distribution at RHIC (0.2 TeV) with 
b = 7 fm. For the ISF eccentricities, both ip\ — ip2 and ip4 — ip2 are mostly around 7r/2 while 



3 The other measurement of 1/% is to project the to the 2nd harmonic event plane, i.e. < v\ cos[n(ip^ — 
1&f)] >. 
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FIG. 6: (Color online). The b dependence of and < cos[n(^ - ip n )] > for RHIC {yfs = 200 
GeV) based on three different models. 



■03 distribute randomly. These seem consistent with Fig. [19] except for the 1st harmonics. 
From — tpi and ipi—ip2 distribution, we expect ipi — ip2 to be around ir/2, but Fig. [8] shows 
ipi — 1P2 to be around 0! This problem is resolved by noticing that the ip{ — ipi distribution 
is broader when ip\ — ip2 is around 7r/2 compared to when ipi — ip2 is around 0. Furthermore 
it implies the conventional picture that ip{ (or ipf) randomly distributes is flawed. 

It is also interesting to separate the contribution to u% from jet spot fluctuation and 
from the matter density fluctuation. Fig. shows the v% due to the full fluctuation (i.e. 
total response) and only due to the jet spot fluctuation (i.e. JS-response) in the central 
collision at RHIC and LHC (2.76 TeV) by using the NTcE model. Here the JS-response is 
computed by removing the fluctuations in the matter density. We can see for the first three 
harmonics, the jet spot contribution is significant. The difference between total response 
and JS-response can be considered as the response due to the matter density fluctuation (i.e. 



11 



0.05 
0.045 

0.04 
0.035 

0.03 
0.025 

0.02 
0.015 

0.01 




LHC(2760GeV) s 

I I I I L 



0.75 

0.7 

0.65 

0.6 

0.55 

0.5 

0.45 

0.4 

0.35 

0.3 

0.25 



b(fm) 



NTcEv 
NTcEc 



10 12 



L c 




o 
o 




10 12 



0.9 




0.024 


0.85 




0.022 
0.02 


0.8 


00 


0.018 
0.016 


0.75 


> > 


0.014 


0.7 


CO 

o 
o 


0.012 
0.01 


0.65 




0.008 
0.006 


0.6 




0.004 



0.009 
0.008 
0.007 
0.006 
0.005 
0.004 
0.003 
0.002 
0.001 




10 12 



FIG. 7: (Color online). The b dependence of and < cos[n(^ - ip n )] > for LHC (y/s = 2760 
GeV) based on three different models. 



MD- response). To compute the difference, we need the information about the sign of the 
responses as shown in Fig. [10] which plots the corresponding i/j^ — ifj n distribution. The plots 
indicate that for the 1st and 3rd harmonics in both experiments, the total and JS-response 
are positive, but in the 2nd harmonics, the total response is positive and the JS-response 
is negative. Combining the information about z/^ with ip^ — ip n distribution, we find that: 
for the 1st harmonics, the total response is dominated by the JS-response while the MD- 
response plays a cancellation role; for the 2nd harmonics, the MD-response is positive and 
the JS-response plays a cancellation role; for the 3rd harmonics, both MD and JS-response 
are positive. For the other higher order harmonics, the JS-response becomes less significant. 
The result is also consistent with our earlier simple estimate jijj]. This shows that the jet 
anisotropy is a probe to both the initial collision density profile at very early stage of the 
AA collision and the mater density profile during the jet traveling through the medium; 
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different u% with the two playing different roles put important constraints on modeling both 
effects. 

Furthermore, is the jet anisotropy sensitive to the composition of the entropy density? In 
Fig. [HI we compare the results for LHC (2.76 TeV) based on two different medium (6 = 
and 7 fm): the "NTcE" calculation is the same as the "NTcE" for LHC (2.76 TeV) in Fig. E 
assuming the two-component profile, S(t ) ~ (1 — 8) X p p /2 + 5 x p c , but the "p p NTcE" 
assumes S(tq) ~ p p . This shows that the first three harmonics are sensitive to different 
entropy densities. Especially, v\ with p p density profile is 60% and 40% bigger than with 
the mixed profile at b = and 7 fm. The ip^ — ip n distributions in the two cases, as we have 
checked, are quie close. Although the conventional procedure uses the multiplicity vs. b to 
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GeV) due to full calculation of NTcE model ("NTcE" curve) and the contribution of the jet spot 
fluctuation in NTcE model ("jet spot NTcE" curve). 

calibrate 5 in the two-component entropy density, our results suggest studying jet anisotropy 
to be another way to constrain the density profile. 

IV. HARD-SOFT CORRELATION 

Motivated by the experimental analysis in Ref. [69f , we study the (un)triggered azimuthal 



54, 69 



correlation between hard and soft hadron at the mid rapidity Following Refs. 
the triggered di-hadron correlation is the pair distribution relative to angle difference, i.e. 
cj) h — s , with (p h constrained in specific region, 1Z. Here TZ is composed of four different 
pieces: <p a — 4>p ^ (p h — ip EP ^ <p a + (ftp (0 a is in the first quadrant) and other three with 



7T + . 



7T- 



In Appendix [HI we derive necessary formula for the computation. 



Similar calculation has been done in Ref. |54j, but we keep track the difference between ip^ 
(ipJ) and V'f ( a l so called as event plane angle ip EP ). To compute the untriggered correlation, 
we can simply set <J) a = <fip = it /A. Because the medium evolution in the transverse plane and 
the subsequent hadronization are not included in our simulation, we simply assume that (1) 
the event plane angle ip EP is the same as (2) v s n = x s n e n for the first four harmonics [§] 4 ; 
(3) the perfect alignment between and ip n . We set Xn=i 2 3 4 5 ~ 0-15, 0.26, 0.21, 0.14, 0.086 



respectively (see e.g. |24j, |27|, |29_|, |30J), which are for the associate hadron with p t around 
2-4 GeV at RHIC (0.2 TeV). 

In Ref. 0], we made a simple estimate of the untriggered di-hadron correlation in the 
central collision at RHIC (0.2 TeV). Here we can improve the previous calculation by sim- 
ulating the realistic density profile and including the ip£ — ip n distribution (ip^ = ip n was 
assumed before). Let's apply the following decomposition: 



C(A0) = 2 x 



^ V n A cos(nA0) 



where the definition of C ( A0) can be found in Appendix [Bl Here we collect the results of 



4 Ref. [30| shows that ~ e„ is valid for the first three harmonics in general but only valid for the 4th in 
the central collision. 
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FIG. 10: (Color online). The distribution of ^ ~ i>n for RHIC (s/s = 200 GeV) and LHC 
(-y/s = 2760 GeV) due to full calculation of NTcE model ("NTcE" curve) and the contribution of 
the jet spot fluctuation in NTcE model ("js NTcE" curve). 



current calculation (previous estimate) 5 : Via = 2.3 (2.8) x 10~ 4 , V 2 a = 3.1 (3.0) x 10~ 4 , 
and V3A = 2.1 (2.3) x 10~ 4 (higher harmonics are at 10% level of the first three). Both 
calculation are based on the same NTcE model. Fig. [12] shows the comparison for C(A0). 
We can see the structures of C(A(j)) are quite close. 

In Fig. [TBI following Ref. [25[, we show our results for V n & vs. trigger direction <p a with 
b = 7 fm at RHIC (0.2 TeV) based on the NTcE model. The 0^ is set as 7r/24, as in 



5 The results of the previous estimate, as shown here, is with parameter £ set as 1.25c. See Ref. [3| for 
details. 
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curve) in Ref. [3j. The two calculations are all based on the NTcE model. £ = 1.25a is a reasonable 
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Ref. [69(| . The even harmonics V2A and V4A are dominated by cos(2A^|)t 2 and cos(4A?/>£ )i 4 
[see Eq. (1B4j) ]. which lead to cos(20 Q ) and cos(40 a ) dependence. The odd harmonics Via 
and V3A are much smaller (in the plot, they are scaled up by a factor of 100). The Via 
has some cj) a dependence because of the mentioned nonuniform distribution of ip{ — ip2 and 
ipi — as shown in Fig. [HI while V3A depends much less weakly on <p a . The author in 
Ref. [25| extracted V n & directly from the unsubtracted correlation data in Ref. 69( 6 . Note 
our Via is positive, but Via of Ref. 25| is negative, which however is consistent with that 



v\ of the associate hadron is negative when p t is below 1 ~ 2 GeV, and positive when p t 
becomes bigger. In addition, the extracted Via can have contribution from the correlation 
due to momentum conservation, although a rapidity gap was required in the extraction. 



6 Although the trigger p t is around 2^4 GeV in the extraction 25| , it was mentioned that for the higher 
Pt trigger, 4^6 GeV, the 4> a and associate pt dependence of V n A are similar as with lower p t trigger. 
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FIG. 13: (Color online). The triggered di-hadron correlation V n A vs. trigger direction <p a with 
b = 7 fm at RHIC (i/i = 200 GeV). (f>p is set as vr/24, following Ref. 69(. The calculation is based 
on the NTcE model. 



Our V 2 a is similar to that in Ref. [25], but V^a's overall sign is opposite to the V4A in 
Ref. [25j, because we approximate cos(4A-?/>£) in Eq. HB4[) by cos[4(-?/>4 — V^)]- It turns out 
although 1^4 — 1^2 is around tt/4 (see Fig. [8]), the hydrodynamic calculation shows for the final 
hadron spectrum ipf — ip2 is around [30j! So we expect in the full calculation with medium 
evolution included, the sign of V4A can be corrected. Because has no p t dependence in 
our calculation, the pt dependence of J4a is the same that of v s w which is consistent with 
the analysis of p t dependence in Ref. [25 . 



V. SUMMARY 

In summary, we have performed a systematic and event-by-event study of the jet az- 
imuthal anisotropy for both Au — Au collision at RHIC (\fs = 0.2 TeV) and Pb — Pb 
collision at LHC (\/s = 2.72 and 5.5 TeV). The Monet-Carlo Glauber model is used to 
generate the initial state density profile that fluctuates strongly from event to event. Three 
different geometrical jet energy loss models, including NTcE, L 2 , and L 3 that represent the 
characteristic path length and density dependence of the near-Tc enhancement, pQCD, and 
AdS/CFT models, are explored and compared at both RHIC and LHC (2.76 TeV). In each 
event, we extract u^, the high p t event plane angles if)£ and the initial participant plane 
angles ip n . For the 2nd harmonics, we see ^2 is strongly correlated with ip 2 i n peripheral 
collisions as predicted by all the models for different experiments. It is reasonable to assume 
a strong correlation between ?/>f and ip2, as demonstrated by the hydrodynamic calculation 
(e.g. |30j). Based on this, we can directly compare the v\ with the measured v 2 of high p t 
hadrons, as carried out in Ref. [4j. There it was found that the NTcE and L 3 models can 
explain the high p t z/ 2 data at RHIC while the NTcE and L 2 models succeeds in describing 
LHC (2.76 TeV) data. One therefore concludes that taking together the geometric data at 
both RHIC and LHC there is only the NTcE model describing all data sets. In addition as 
discussed in Appendix D in detail, the NTcE model naturally explains the less opaqueness 
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of the medium created at LHC as compared with RHIC, i.e. a reduction of jet-medium 
interaction at ~ 30% level as implied by data. We also analyze the correlation between ip^ 
and ip n for other harmonics. For the 1st and 3rd harmonic, < cos[n(^ — ipn)} > are around 
0.6 ~ 0.8 except in very peripheral collision, but for the 4th and higher harmonics, even the 
sign can change although the specific transition impact parameter depends on models and 
experiments. This seems consistent with the LHC (2.76 TeV) data shown in Ref. 0] which 
right now are not converging but does indicate negative responses for higher harmonics. We 
further clarify that in the middle-centrality collisions (e.g. 6 = 7 fm) at RHIC, although 

and ^3 does have a random distribution, the 1st harmonics turns out to be tricky: ipi 
has slight preference in the out of plane direction but ip{ prefers to be along the in-plane 
direction. For the 4th harmonics at 6 = 7 fm, both -04 and fluctuate around 7r/4 relative 
to either participant plan ip 2 or the reaction plane. It is also interesting to separate the 
contributions to the jet anisotropy. In the central collision at RHIC, we have analyzed the 
effect of the initial jet spot fluctuation. The difference between the JS-response and the 
full response should be due to the matter density (shape) fluctuation. We see in the first 
three harmonics, their roles are different, which points out the importance to have a coher- 
ent picture for all these harmonics in a valid model. We also test the sensitivity of the jet 
anisotropy to the composition of the matter density in the Glauber model, which shows at 
LHC (2.76 TeV), the first three harmonics are affected significantly by changing the entropy 
density from a combination of participant and collision density to the participant density 
alone. Although the mixing of the two in the matter density is normally calibrated to the 
observed multiplicity vs. 6, our test does show another interesting way to calibrate this. 

Based on the information about jet anisotropy and assuming that the response of the 
collective flow to ISF is linear, we proceed to discuss both untriggered and triggered hard-soft 
correlation. For the central collision at RHIC, we improved our calculation of the untriggered 
correlation in Ref. [i| and confirmed the results there: a strong peak develops on the near 
side ("hard ridge") while a double-hump structure shows up on the away side. This could 
be a possible explanation for the experimentally observed "hard ridge", but the away side 
structure can be more complicated because of the other sources of correlations (e.g. effects of 
global momentum conservation and cluster correlation 70, 7l|). which requires more detailed 
study. By using our re-derived formula for the triggered correlation, we have studied the 
Fourier component of the correlation, V n ^ . We find that the trigger-angle dependence of V n ^ 



is consistent with the extracted information in Ref. [25(. As for the associate p t dependence 
of the hard-soft correlation, we expect it to be the same as that of z^(pt), because of the 
weak p t dependence of z/^ as shown in data and implemented in our simulation. This also 
agrees with the findings in Ref. [25] . However a decisive comparison with data again requires 



a good understanding of other sources of correlations. 

We end with discussions on a number of issues to be improved in future studies. A 
detailed study of the hard-soft correlation naturally requires including a realistic medium 
evolution in the transverse plane through hydrodynamics j72|, 73]. The linear response is a 



crude approximation for the soft dynamics. Moreover, such transverse expansion dynamics 
can also play a role in the jet quenching, in particular for higher harmonics. In this simula- 
tion, the transverse shape of the medium is frozen (matter density decreases as 1/r because 
of the longitudinal expansion), but in reality the eccentricities of the matter will reduce with 
increasing time because of the pressure gradient driven expansion (and such "self-quenching" 
may be a significant factor for higher harmonics) [t3] - Although qualitatively, the major- 
ity of jets may not experience this change very much, because the dramatic shape change 
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happens in relatively later stage, it is certainly interesting to explore this effect especially 
for harmonics other than the second harmonic, which are mostly driven by the shape fluc- 
tuations. Integration of our current jet azimuthal anisotropy study with hydrodynamics is 
being pursued. Furthermore, the p t dependence of jet quenching is not implemented in the 
present geometric models but efforts are underway to include that. Another important un- 
certainty is related with possible pre-equilibrium energy loss, which hopefully will be better 



understood and estimated with improved descriptions for the pre-equilibrium evolution [75 



We are focusing on the geometrical information about jet anisotropy, but it is important to 
have a coherent understanding in terms of different variables, which requires a careful study 
of the hard parton energy loss for example in the near-Tc enhancement picture. Finally, we 
hope we have made a case here for using the hard probe as a new and sensitive tool for quan- 
tifying the initial fluctuations in heavy ion collisions and for discriminating different models 
for the initial conditions. Future studies along this direction will be certainly extended to 



hard probe of initial conditions generated from e.g. CGC-motivated models |76l . 177] and 
their comparison with MC Glauber model. 
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Appendix A: The structure of the simulation code 

Fig, mi summarizes the overall structure of the simulation in each event. It starts with the 
sampling of nucleon positions in the two colliding nuclei. Based on the Glauber model, we 
can identify the binary collision pairs, which leads to the participant density and collision 
density by using proper smearing procedure. This provides the basic information for all 
the following calculation modules. In the jet quenching module, we calculate the azimuthal 
distribution of the Raa as described in the text and get v% and ipn- The 2nd module analyzes 
the initial entropy density, which gives e n and ip n . In principle other modules can be included 
in the code, for example the EM field calculation. Then the results from different modules 
are grouped together and recorded in the output file. To get an ensemble of events, the 
collision geometry, e.g. impact parameter for AA collision, needs to be sampled, which is 
not shown here 7 . After the initialization of the collision geometry, the whole calculation 
shown in Fig. Ocan be carried out. 



7 In U — U collisions, the collision geometry is more complicated, requiring the sampling of impact parameter 
and the relative orientation of the two nuclei. 
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FIG. 14: The structure of the simulation code. 



Appendix B: The formula for the triggered di-hadron correlation 

In this section we derive the formula for computing the triggered di-hadron correlation 
mentioned in Sec. IIV} based on the following soft and hard hadron azimuthal distribution: 

dN s 

— - ~ 1 + 2 cos[n(0 s - ^ s n )} , (Bl) 

n 

a — ~ 1 + 2 £ v h m cos[m(0* - rj] ■ (B2) 

m 

According to the definition in Sec. IIVt let's calculate the following: 

/dN s dN h 
, d(t)h ~d^l^ 5{(t>8 ~ ^ ~ A(t>) ~ 1 + C(A<t)) ^ 1 + 2 ^ C ° S(nA0) ' (B3) 

Here 1Z is composed of four different pieces: 4> a — (ftp ^ 4> h — ip EP ^ 4> a + 4>i3 (4>a is in the 
first quadrant) and other three with <j> a — > —<f) a , n + <p a , 7T — <p a . We keep track the difference 
between (ip£) and ?/> EP , which is different from Ref. 54]| . By using Eqs. f)B2j) and fIBlj) . 
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If 



In the above expression, Aip£ = ip£ — ip ] Aip% = ip^ — ip 

we assume all the reference angles ip^ and i\) 3 n are exactly correlated with ip EF , the above 
formula would become the one in Ref. 0- To compare with the data, we need to average 
C(A(p) over many events, i.e. < C(A<f)) >. In principle, C(A0) should also have sin(nA0) 
components, but their amplitudes should be very small with many events averaged. [We 
will simply use C(A0) as < C(A0) > from now on.] 



Appendix C: Discussion on the near-T c enhancement and the trace anomaly 



We discuss here the possible connection between the near-T c enhancement of jet-medium 
interaction and the QCD trace anomaly (e — 3p)/T 4 (with e,p, T the energy density, pressure 
and temperature) which as an interaction measure also develops a strong near-T c peak as 
shown by lattice QCD simulations [' 



7$. 



One may naturally ask whether these two 
non-perturbative effects near T c could share a common underlying picture. We notice 



that there have been attempts [801 to develop non-conformal holographic models based 
on gauge/gravity duality that can mimic the QCD trace anomaly well but face tensions 
in describing jet quenching phenomenon. In this appendix, we will show that the near-T c 
enhancement and the trace anomaly may be consistently described together by a plasma of 
magnetic monopoles based on the "magnetic scenario" in 0, [13, H S3 



We consider QCD near-T c plasma as an ensemble of thermal magnetic monopoles with 
certain density n m {T) and effective mass M m {T). Let us first compute the contribution to 
trace anomaly from this magnetic component. As a crude approximation let us first neglect 
the potential energy and express the trace anomaly through the density, i.e. 



e — 3p 

rp4 



Xm 



(CI) 
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with the dimensionless coefficient xm defined as: 

x 2 dx M 2 



/T>2 \ ° exp(\Jx 2 +M 2 \-l \fx 2 +M 2 

Xm = ( r=„ =) = J Ax (C2) 

■\Jx 2 +M 2 



:x' + m 2 / Jo — 

exp 

In the above the thermal monopole mass M = M m /T is about 1.5 ~ 2 when close to T c 



83l . |84| . which corresponds to xm ~ 0.7 ~ 1.1 . We therefore have roughly the relation 



3 " ~ (0-7 ~ 1.1) x ^ (G3) 



Now the monopole density close to T c , according to lattice results and model studies [1711811 



84j . is about n m /T 3 w 2 ~ 3, so we get an estimate of the monopole contribution to trace 
anomaly near T c to be about 

^ « 1-4 ~ 3.3 (C4) 

Note that the peak value of trace anomaly found from lattice simulations is about 

so the above contribution is well consistent with lattice data and 



actually is the dominant part. The mild difference may be made up from the interaction 



among monopoles: according to the analysis in [16|, [83[ the interaction potential energy 
is comparable with the kinetic so could contribute similar amount to the trace anomaly. 
We also note that the rapid decrease of trace anomaly with increasing temperature can be 
understood from the rapidly dropping density n m /T 3 for T > T c . 

Let us now turn to the discussion of the jet-medium interaction within the same magnetic 
scenario. With similar picture to the GLV model for jet energy loss \^ while replacing the 
Coulomb electric scattering centers by the magnetic ones in the near-T c plasma, we expect 
roughly the amount of medium kick received by a penetrating electric jet to scale as 

~ a E (T)a M (T)n m (T) (C5) 

While the electric coupling ag and magnetic coupling a M both depend on temperature scale, 
their product has to be unity as per the famous Dirac condition, i.e. «e(T) omI? 1 ) = 1. 
Therefore, we see that if we normalize the medium kick by the entropy density at the same 
temperature, we will get 

a E (T)a M (T)n m (T) n m /T 3 

K W) = ^jT r (C6) 

Again the monopole density n m /T 3 strongly peaks at T c (T^. 



83L [84j while s/T 3 relatively 



slowly decreases from high temperature to low temperature [78|, [79|. As a result the normal- 
ized jet-medium interaction strength k develops a strong near-T c enhancement. To make a 
more quantitative statement, let us compare the k at two temperatures: from 1.5T C to T c 
the monopole density increases by a factor of ~ 3 while the entropy density s/T 3 decreases 
by a factor of ~ 2, and therefore the k increases by a factor of ~ 6 from 1.5T C down to T c 
which is well in accord with the assumed enhancement in NTcE model. 

To sum up the above qualitative discussions and estimates, we see that the presence 
of an emergent monopole plasma near-T c with the density and thermal mass indicated by 
lattice and model studies, can consistently describe the strong near-T c peaks in both the 
trace anomaly and the jet-medium interaction. 
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Appendix D: An estimate of average opaqueness from RHIC to LHC 



In this part we make a simple estimate, in the NTcE model, of the evolution of the 
average opaqueness of the fireball from RHIC to LHC. The hot medium created in heavy 
ion collisions is neither homogeneous nor static, and therefore one can only talk about the 
opaqueness of the fireball created in the collision on an average sense, i.e. averaging both 
over the spatial distribution and the time evolution. For simplicity let us consider the 
perfectly central collision with 6 = 0. We will use this case as an indicator of the opaqueness 
evolution with collision beam energies. The average opaqueness seen by a jet going through 
the medium along a particular path P can be estimated as follows: 

J p n[s(l)]s(l)ldl 

< K[S{1)] >P= J p s(l)ldl (D1) 

Note that even for the same geometry and same jet path in the fireball created at different 
collision energies, the s(l) will be very different and the convolution with re[s(Z)] is rather 
nontrivial due to the peculiar near-T c enhancement form in Eq.(jSJ). This is of course further 
subject to averaging over paths, i.e. over all initial jet spots weighed by collision density as 
well as all possible jet directions. We will use a simple optical Glauber model (with the same 
parameter sets as in the previous MC Glauber simulations) to evaluate this average. Note 
that the longitudinal boost-invariant expansion effect is also included as in the previous MC, 
i.e. the density s(l) on a specific point of the path is really a function of both I and time t. 

To give a quick idea, we can simply consider a jet initiated right from the center r± = 0: 
with this simplest path we can easily obtain < k >rhic : < K >lhc ~ 1 : 0.75. Clearly we 
see a shift toward less opaqueness as a consequence of two factors: the near-T c structure of 
k(s) and the shift to higher density matter at LHC. In contrast, all models with a constant 
k will see no change in the opaqueness from RHIC to LHC. A full evaluation averaging over 
initial jet spots and paths with fluctuating initial condition gives the following evolution of 
average opaqueness: 

< « >rhic:< k >lhc~ 1 : 0.72 (D2) 

We note this is not only qualitatively but quantitatively in agreement with the earlier 
analysis in Q which found that in order to describe the LHC Raa data the jet-medium 
coupling parameter k (which is a constant in the model of Q) has to be reduced by a factor 
~ 30% from the RHIC value in the same model. What we want to emphasize here is that 
such a reduction of (average) jet-medium interaction implied by data from RHIC to LHC is 
naturally born out from the strong near-T c enhancement of jet-medium interaction! 
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FIG. 15: (Color online). The distribution of ip n with 6 = 0, 7 fm at LHC (y/s = 2760, 5500 GeV). 
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FIG. 21: (Color online). The distribution of 
based on three different models. 



i) n for LHC {yfl = 5500 GeV) at b = 0, 7 fm, 
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